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1. List of Symbols 
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M 
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c 

c 

I 

c 

r 
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E 
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G 

g 

h 


i/ j/k 


1 
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M" 

M 

M 

y 
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Integration points in the negative Z direction 
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n 

N] 

N, 

I 

N: 

N 

I 
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NW(1) 

NW(2) 


P' 


P 


t 


Q 
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R 


R 


I 


R 


V 


^3 
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1 


S 


S' 


Integration points in the positive Z direction 

Unit vector normal to wing surface 

Upstream integration points 
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chordwise direction 
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chordwise direction 
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Perturbed pressure 
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2 
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Spanwise integration region away from the wing 
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2. Introduction 


When the freestream is nonuniform, the assumption of 
flow irrotationality becomes inapplicable so that the 
conventional potential flow theory must be revised (ref. 1). 
Practical examples of nonuniform freestream include the wing 
behind a canard wake, a wing or tail situated in the propeller 
slipstream and airplanes in a wind shear. 

The literature on the subject is not extensive. Von 
Karman and Tsien (ref. 2) developed a lifting line theory for 
an incompressible nonuniform stream in 1945. Homentcovschi 
and Barsony (ref. 3) formulated a lifting-surface integral 
equation for nonuniform incompressible flow with a general 
stream velocity profile. Hanin and Barsony-Nagy (ref. 4) 
developed a slender wing theory for a subsonic or low- 
supersonic nonuniform stream. Later, a lifting line theory 

(ref. 5) for wings in a subsonic nonuni^orm stream was 
presented. Recently, a lifting surface theory for nonuniform 
supersonic parallel stream (ref. 6) was also developed by them. 
More recently, experimental investigation of wind shear effect 
on airfoil aerodynamic characteristics has also been conducted 
(ref. 7). K. Gersten and D. Gluck (ref. 8) investigated the 
effect of wing wake on tail characteristics theoretically and 
experimentally. 

In the present investigation, the small-disturbance steady 
subsonic flow equation for rotational flow is to be solved. 

The resulting equation is a partial differential equation 
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with non-constant coefficients. It is transformed into an 


integral equation through Green's theorem. It is assumed 
that (1) the Mach number profile M(y,z) or the velocity 
profile has a nonzero value M(y,o) on the wing plane; (2) 
there exists a finite second order derivative M"(y,o) on the 
wing plane, and (3) M' (y,z) is integrable across the stream. 
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3. Theoretical Development 

3.1 Mathematical Formulation 

The partial differential equation for small-disturbance 

subsonic steady rotational flow may be written as (ref. 1) 

„ M (y,z) M (y,z) 

[1-Mco (y,Z)]P' -2 --7 -T~ P'-2 rj-j — — T— P' + P' + P' =0 (1) 

XX M(y,z) y M(y,z) z yy zz 

where M(y,z) is the undisturbed free stream Mach number and P' 
is the perturbation pressure. Here, the flow is assumed to be 
steady, invicid, and compressible past a thin wing at a small 
angle of attack. 

The wing lies in the x-y plane with the positive x axis 
being streamwise along the wing center line. The origin of 
the rectangular coordinate system is assumed to be at the 


wing moment reference point. Let 


X = X , 


= 3 y, z = 6 z, 3 ^ = 1 -m" 


( 2 ) 


After the coordinate transformation, equation (1) becomes 


P' , , + P' . , 

X ' X ' y y 


+ P' . , 

z ' z ' 


M 


= 2 


zl 


P',+2 

y 


^z ' 1 2 —2 

P:. + i-(M^-M^)P;,, 

3^ 


M 'y' “ M "z' „2 “ '"x'x' (3) 

where M = constant and is assumed to be a reference uniform 

Mach number. To use Green's formula: 

^ ^ f^|^)ds (4) 

V g dn c^n 


let 


n = p ' 

^ = i = 


a/ (x-?) ^+3"(y-n) ^+3^" (z-c) ^ 

Consider an arbitrary body with boundary surface S in a 

compressible flow. Outside the body, the whole flow field 

has the volume V. Now, let P^(x,y,z) be any point in V. 


(5) 

( 6 ) 
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Then i|; is a solution of Laplace's equation: 


9^ijj « 

^ ^ 9^ - ° 


( 7 ) 


However, the basic solution becomes infinite at so that 
Green's reciprocal formula can not be applied. To avoid 
this difficulty, surround with a small sphere S' with 
radius e. Let V be the volume which is obtained by excluding 
from V the interior of S'. V is bounded by S and S'. 

Applying Green's reciprocal formula to equation (3), it is 
obtained that (see Appendix A) 

P(x' ,y' ,z') 

i'lg Ik'- Ik (^)ia5an-2J'fp'dcdn 


+2 


; ; ; M P'd5dnd^-h2^^/_9ji p'd^dndc 


V 

_ ;;; l 


9n 


V 


(M^-M P'^^d^dndc} 


( 8 ) 


The boundary condition (the flow tangency condition) is given 
by (ref. 1) 


/ \ 

w(x,y , z) = c - a 

V(y,z) 9x 


(9) 


Following the linear theory, the boundary condition is satisfied 
on Z = 0 plane. 

From reference 1, the linearized momentum equation in Z 
component is 


9w 


9p 


Pco(Y/Z)V^(y,z) = 0 


( 10 ) 


It follows that the downwash in equation (9) can be written as 


w(x,y,z) = 


-1 


9^ (x.,y,z)dx. 

Poo(y/2)V^(y/Z) 


( 11 ) 


where the lower limit is chosen in such a way that w(-°°,y,z)-0. 
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Using equations (8) - (11) , after lengthy manipulation, the 
following final integral equations and boundary conditions 
for two-dimensional and three-dimensional flows can be 
obtained (see Appendix A and Appendix B) . 

(a) Two-Dimensional Flow 

The integral equation is given by 


C^(x,Z^) =-^ j Z„ACp 


2tt 0 [ (x-5) ^ + 6"Z^M 


dC 


- 1 S^(0) 


2tt6 M 


r' Jin[ 3' ■hC(C)dC 

0 1- (x-a 


M/ 


+ 1 


2tt3 


{ 


23^M (Z,,-C) 


[ (x-C) ^+3^ {\-0 ] 


+ B(c) £n f 3" ^]}c (C,rjd?dc 

L(x-a "+3" (Z^-C) "i 

+ 1 ff Sg(g,c) (x-g) (M"-M^) dgdc 

2 ^ [ (x-U ^+3^ (Z^-c) 


The boundary condition is 

\v (x,0) =_^c - a 
9x 

= -Bf^ ACp(g) dg+ • B(C) (tan~ ^ x-C +tt/ 2 ) 

4it ° x-C 27 t I C I I -3C I 

Mr 

+ (x-5) 1 C_(5,;)d5d5 

+ ; j (mMc)-m^)-c» dgdc 

4tt 3 [ (x-C) ^+3^C^] 

9 Mr 3 

Where B(^) = )/ ^ = Mach number, = y^(M), and 

Cpf = (?/?)• is the pressure coefficient. 

S o s P P 


(12) 


(13) 


in 



(b) Three-Dimensional Flow 


The integral equation is given by 


Cp(x,y,Zv) 


4tt 


; I 


$ 2 Zv-ACp(?,n) 


- [ (x-^) ^ (y-n) ^ + 3^z ^ ] 


3/2 


d?dn 


’w 


j I 


Mr 

^(n, 0 ) ACp(c,n) 


[ (x-C) " + 6My-n) " + 1/2 


d?dn 




, Mr 

^ M • (Zv"C) 


V 


, [ (x-c) ^ + 3My-n) ^ + 3^ (z^-c) 


jj2 * Cp(C/n/C)d^dndc 




B(n,^:) 


V 


. [ (X-?) 2 + 3My-n) ^ + 3 Mz^-c) 


jj2 ‘ Cp(C/n/C)dCdndc 


1 //; 


2 Mxi 

3 (h/C) * (y-n) 


2T7 


V 


, [ (x-u " + 3My-n) ' + 3"(Zv-0 M 


jj 2 * Cp(C,n,^)d^dnd^ 


+ III 

2tt 


c (n^c) 


[ (x-5) 2 + BMy-n) ^ + 3 MZv-?) 


2^^ 2 * d^dnd^ 


(x-^) ♦ (M^-M ) 


I I I 


y, [ (x-C) ^ + 3My-n) ^ + 3 MZv-C) M 


■jy^ • Cp^(C,n,c)dCdnd? 


( 14 ) 
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The boundary condition is 
W(x,y,z) = 


9x 


- a 


8it 


2r7 2 


= M AC„(g,nM ^ 'y-"* ^ (1 + 


X - ; 




[3^ (y-n) ^+3^z^ ] ^ [ (x-^) ^+3^ (y-n) ^+3^z^] 

3^zMx-g) 


1/2 


) 


M; 


. Bi ,, 

s„ 3^ (y-n) ^+3^z^ 


[3My-n) 2+3 ^zM [ (x-C) ^ + 3^ (y-n) ^+3^zM 

x-i 


3/2 


"W 


?! 


1 + 


[ (x-^) ^ + 3" (y-n) " + 3^z"] 


1/2 


d^dn 

d^dn 


- ^ ''' Cp(C,n,c) 


V 


C 1 + 


M_ 

-^[ (y-n) '-(Z-?) 2] -B(z-0 [(y-n) '+(z-c) 2] 

[ (y-n) " + (z-c) 2] 2 
T72 ) 


X - g 

[ (X-?) 2 + 02 (y_^) 2^g2 (2-^) 2] 


Mr 

3^ (z-c ) 2 (x-^) 


[ (y-n) ^+(Z-C) "] [ (x-a '43' (y-n) '+3' (z-c) '1^^^ 


d^dnd^ 


!u c^{^,n,o 


4tt ' " ''P 

V 


-C(Z-C) [ (y-n) '+(Z-U^ ^ ^ (y-n) (z-c) 


1 + 


[ (y-n) ' + (Z-c) '] ' 
X - C \ 


[ (X-?) 2 + 32 (y_T^) 2 + 32 ( 2 _ 2 ) 2 ] 1/2 


Mn 

3' (x-^) (y-n) (Z-?) 


[ (y-n) '+(Z-C) '] [ (x-C) '+3' (y-n) '+3' (z-c) j 


d?dndc 


^ (5.n,U 

v' [ (x-C) '+3' (y-n) '+3' (z-c) '] 


372 


(15) 


where 


9 Mr 

B(n,0 - - (-^) 


9 

and C(n,^) =-^(-^) 
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3.2 The assumed functional form for Cp 
3.2.1 Two-Dimensional Flow 

The unknown pressure coefficient Cp is assumed to have 
a functional form which can be derived by retaining only the 
first term on the right hand side (R.H.S.) of equation (12). 
It follows that 


Cp ( ^ , i; ) 
Cp^ (C,C) 


K 


2 

j=i 


? 


3? 


- tan 


-1 ^ 


3C 


] 


K 

V A- 
3 I? 

j=l 


[ 


Be 


Be 


B"e'+(C2^-?) ' 


.2^2 


e^+(?,j-a 


■] 


(16) 

(17) 


where £, . and £ . are chordwise control point locations given 

by a cosine law distribution as illustrated in Figure 1. 


3.2.2 Three-Dimensional Flow 

In this case, the assumed functional form for Cp is 
derived by retaining the first term on the R.H.S. of equation 
(14) with etc., defined in Figure 2. 


Cp ( 5 / n / e ) 


Ny N2 

j. Y Y A 
2tt Z. ij ? 

j=l i=l 






i— _ (tan-^- 

c I. I V 


(n-n^.) U-C^i 


==) 




elJ(C-? 2 i) "+B" (n-n^j) '+B"e 

r -1 I 

' 1 elv/(?-C'i) ^+BMn-n,j) '+B"e' 

^ii“^ ^ _1 (n-ri2j) 

tan . 


+B^e^) 


5I J(?-?;^) 2+BMn-Ti2j) '+BH 

-iiAli- (tan-^ ' ) 

(18) 
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Cp 


Nz Ni f 

= -^ Y i 

n~n • / 
n / 

1 


\ 

2tt Z< ^ ijC) 
j=l i=l L 

_l+Di^ ' 

>/("c-C^i) '+Bi 

[(C-5^i)"+Bi] 

3/2/ 



1 


\ 


l+Dz" ^ 


[ (5-^2i) "+Bzl 

3/2 ] 


ri-h • / 

n ( 

1 

(C-^i)^ 

1 


. l+Ds" ' 

y("^-c;i) "+Bi 

i(5-e;.)"+B2) 

3/2 ) 


riii 1 

1 

'E-E>i)^ 

ll 


1+D4^ ' 

7( E”E ^ ^+B2 

.) "+Bz] 

3/2 i 

where 




(19) 

1 

1 

2i'lE-q 

il 




D 2 — 


Da = 


D4 = 


Icl/IF 

•5U> 

" + B" (n-ri2j) ' 

+ B^?^ 


(n-n 



Ul7(T 


" + BMn-n^j) " 

+ B'c" 


(n-n 

.i'lC-EJil 


UI>/(? 


" + BMn-n^j) ' 

+ B"?' 


(n-n 




K 1 y ^ (n-rij j) ^ + 32 ^- 


Bi = B' (n-n^j) ^ + B'c' 
Bz = B^ ^ 
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4. Numerical Calculation 


4.1 Solution procedures 

For airfoil problems, equations (12) and (13) are to be 
solved simultaneously. Similarly, for three-dimensional 
cases, equations (14) and (15) are to be solved. 

To simplify the equations, those double integral terms 
in equation (15) are reduced to finite sums following the 
quasi-vortex lattice method (ref. 9) . To satisfy the wing 
boundary condition, the continuous vortex distribution over 
the wing is replaced by a quasi-continuous one, being 
continuous chordwise but stepwise constant in the spanwise 
direction. Thus, the wing surface can be divided into a number 
of vortex strips with the associated trailing vortices (Figure 
3) . In any strip, consider a vortex element y dx with an 
arbitrary direction Z (Figure 4) . The integrals are then 
reduced to finite sums through the mid-point trapezoidal rule 
(ref. 9 and Appendix B) . Similarly, those double integral 
terms in equations (14) can be simplified in the same way 
(see Appendix B) . 

4.2 Singularities and integration regions 

Before proceeding with the numerical integration of 
integrals, singularities of integrands must be examined. 

From the definition of 4* in equation (6) , it is expected 
that all integrands will have a singularity at the observation 
point x~F,, y=rif z=C. To account for this singularity 
properly, rectangular coordinates are transformed into polar 
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coordinates by using a cosine law relation (i.e., half circle 
relation) with separated control and vortex points (ref. 9). 
Outside the wing, a quarter-circle transformation is used 
which creates smaller mesh-size near the airfoil or wing. 

The mesh size is small enough to represent better the rapid 
change in the pressure coefficient. The detailed integration 
regions and coordinates transformation are presented in Appendix 
C. 

4.3 Numerical Convergence 

The convergence of numerical integration is checked in 
two ways. One is by comparing integrated results of each 
term in equations (12) , (13) , (14) , and (15) , and the other 

is by comparing calculated aerodynamic characteristics. 

4.3.1 Two-dimensional Flow 

The integration schemes for an airfoil are illustrated in 
Figure 5. Calculated ACp's for different integration schemes 
are presented in Table 1. It is seen from Table 1 that 
satisfying equation (12) at 100 points of through the 
least square method produces results for AC which can be 

ir 

obtained by satisfying equation (12) only at = 0.05 (Case 
No. 8 in Table 1) . Note that the airfoil chord length is 
taken to be unity. Therefore, to save computing time, equation 
(12) will be satisfied only at = 0.05 from now on. For 
a wing, this is revised to be Z^ = O.05 Cq, where C^ is the 
mean geometric chord. 

It should be noted that the so-called least square method 
for satisfying equation (12) is based on the following concept. 



By choosing N points of (x, Z^) in the x - z plane, equation 
(12) can be integrated to result in N values of differences 
between both sides of equation (12) . Let these values be 
denoted by F, . Let Fu = ^ b . Ac + ZC. A. (20) 

^ i ^ Pi i ^ 


where K = 1,2. . . .N . 

2 

If the sum of Fj^ is differentiated with respect to A ^ , and the 
results are set to zero, it is obtained that 


v—1 9 A . 


0, j = 1, J 


( 21 ) 


This results in a set of simultaneous homogeneous equations for 
Cp's and A's: 


zc 

k 


jk^ik 



+ 




A. 


X 


0 


( 22 ) 


4.3.2 Three-dimensional flow 

The spanwise integration regions are illustrated in Figure 
6. Using an elliptic wing of AR = 10.91, convergence of 
numerical integration of integrals in equation (15) has also 
been investigated. Note that all integrations are performed 
through the midpoint trapezoidal rule after the coordinate 
transformation described in Section 4.2. From this study, it 
is determined that the following values for integration para- 
meters are appropriate: N^=6, N =0, N„=5, M=5. M^=2, 

i P 3 y ’ 1 * 

=2.5 X C^, = 0.8 X C^, Yg = 0.5 X b/2, and = 

3.0 X b/2 . 
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4.4 Freestream profiles 

In the present investigation, the freestream velocity 
profiles are of three types only — jet, wake, and linear shear 
profiles. These are illustrated in Figures 7 and 8. The 
jet or wake stream profiles are described by 

M(Z) = + (Mq - M^) exp (-(z-z^)^/H^) (23) 

On the other hand, the linear shear profile is given by 
M = for z >. H 

= 1/2 (M^+M 2 ) + 1/2 (M^-M 2 ) (z-z^)/H for -H£z^H 

= for z£H (24) 
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5. Numerical Results and Discussions 

In this section, some numerical results by the present 
method will be presented. Comparison with theoretical 
results from references 5 and 10 will be made first. 

5.1 Two-dimensional Joukowski Airfoil in a Jet Stream 
Results of the present thin airfoil theory for a 2-D 

Joukowski airfoil as shown in Table 1 are compared with 
those obtained by a finite difference method (reference 10) 
in Figure 9. It is seen that the present predicted pressure 
peak is more aft and the peak magnitude is slightly less than 
those given by reference 11. Thickness effect may be 
responsible for the discrepancy between these two methods. 

5.2 Elliptic Wing in the Jet and Wake 

Elliptic wings of various aspect ratios have been 

extensively investigated in reference 5 by the lifting-line 

approach. The present results for a wing of AR ^ 

compared with those from reference 5 in Figure 10, in jet and 

wake flow. It is seen from Figure 10 that the agreement between 

two theories is good except when is large in the wake 

flow. From the general linearized partial differential 

equation for the rotational flow (eqs. 3 and 23), it is 

1/2 

known that not only (1-M^) AR and are important 

similarity parameters, but also is M^. The theory of reference 
5 does not show the dependence of results on independently. 
The ordinate in Figure 10 is C^/C- , where C_ is the lift 

J_i Li^ J-i 

O 

coefficient in nonuniform flow based on the local velocity 
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V (o) and density p (o) . is the lift coefficient of the 

^o 

same wing at the same angle of attack a in uniform stream. 

The results also show that the effect of a nonuniform 
stream on the lift and induced drag becomes stronger with 
increasing ratio of maximum to minimum Mach numbers of the 
stream. 

5.3 Effect of Wing Wake on Tail Characteristics 

As shown in Fig. 11, a tail surface is frequently 
situated in a non-uniform flow field caused by the wake of 
a wing. The experimental results for a wing-tail configura- 
tion shown in Fig. 12 were obtained by W. Siegler (ref. 13). 

In the present calculation, the wing-section drag coefficient 
is needed. Its assumed values are shown in Fig. 13 by extra- 
polating available data for NACA 0012 airfoil which was used 
in the experiment. The relation between wing and tail angles 
of attack is shown in Fig. 14 and the wake location in Fig. 15. 
(See Appendix D for a detailed calculation of the wing wake). 
The results by the present method for this configuration are 
shown in Fig. 16 together with experimental data. 

In Fig. 16(a), it is seen that with the wing at low 
angles of attack, the wing wake is far below the tail and 
has no effect on the tail. With increasing angles of 
attack, the lift of tail increases until a certain angle 
of attack (24°) is reached. Then the tail enters the wake 
center of the wing and the tail lift decreases as the 
angle of attack is increased further. 

Fig. 16(b) shows that the predicted pitching moment 
exhibits the same trend as the experimental data from reference 13. 
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However, the predicted magnitude is too negative, even in the uniform 
flow. At a certain angle of attack, the tail enters the wing wake, and the 
nose-down contribution of the tail to the total moment is reduced, resulting 
in unstable pitching moment characteristics. 

5.4 Rectangular Wing of AR = 3.3 in the Wake 

Theoretical longitudinal aerodynamic characteristics of 
a rectangular wing of AR = 3.3 in the wake stream are 
shown in Figure 17. In Figure 17(a), it shows that as the wake 
center is far away below the wing, the effect of reduced 
dynamic pressure is small. Due to the positive velocity 
gradient, the vorticity is in the same sense as the circulation 
around the wing. As a result, there is a gain of lift. When 
the wake center moves closer to the wing, the effect of 
reduced dynamic pressure becomes more important and the lift 
is reduced. After the wake center moves up and away from the 
wing plane, the negative velocity gradient contributes a loss 
in lift compared to the lift gain when the wake center is at 
the same distance below the wing . As the wake center moves 
further away, the lift will get close to results of the 
uniform flow. 

In Figure 17 (b) , it is seen that as the wake center stays 
below the wing, the positive velocity gradient makes a significant 
contribution to the leading edge thrust. At Z^=-0.7071, the 
induced drag is calculated to be negative due to large positive 
velocity gradient. Whether this is possible in reality requires 
further study. On the contrary, the leading edge thrust is re- 
duced when the wake center is above the wing. 

In Figure 17(c), the aerodynamic center is shown to shift 
forward (i.e., more negative when the wake center is 
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below the wing and moves backward when the wake center is 
above the wing. Again, the velocity gradient effect is res- 
ponsible for this difference. 

Figure 17(d) indicates that the pitching moment becomes 
less negative as the wing moves into the wake. This trend is 
consistent with the results presented in reference 8. Within 
a certain range of angles of attack, as the wake moves toward 
the wing, 9C^/9a becomes positive, contributing to pitch insta- 
bility of the airplane. 

5.5 Rectangular Wing (AR = 7.2) in the Linear Sheared Flow 

Theoretical longitudinal aerodynamic characteristics 
of a rectangular wing of AR =7.2in the linear sheared flow 
are shown in Figure 18. 

In Figure 18 (a ) , it is seen that the lift is decreased 

due to the local dynamic pressure effect. However, the 

positive velocity gradient makes a positive contribution to 

lift as mentioned earlier. With Z, = 0.0 and M =0.15, the 

d o 

lift is slightly greater than those in uniform flow (Moo = 0.15) 
due to this positive velocity gradient effect. For the sheared 
flow, since the second derivative of the Mach number profile, 
M"(?)/ is zero, the change of lift and induced drag coefficients 
is much lower than those in the wake or jet. 

In Figure 18(b), since the velocity gradient effect is 
almost the same for each case, the difference in lift-drag 
ratio is mostly due to the local dynamic pressure effect. 
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As shown in Figure 18(c), the aerodynamic center is 
not much changed in the linear sheared flow from that in the 
uniform flow. This is because nonuniform flow effect makes 
the same contribution to both the lift and pitching moment. 

5.6 Plane Delta Wing (AR = 1.4559, A= 70°) in the Jet 

It is assumed that the suction analogy of Polhamus (ref. 

12) is still applicable in a nonuniform flow. The vortex lift 
calculation through the method of suction analogy described 
in reference 11 is applied to a wing in the nonuniform free 
stream. The longitudinal aerodynamic characteristics of a plane 
delta wing of 70-degree sweep by the attached flow theory and 
by the vortex lift theory are shown in Figure 19. It is seen 
from Figure 19(a) the vortex lift increase is slightly larger 
in the jet stream than that in the uniform stream. Figure 
19(b) shows that the lift-drag ratio is increased in the jet 
either in the attached flow or in the vortex flow as compared 
with that in the uniform flow. This is because the wing is at 
the jet center, so that the local dynamic pressure effect 
plays the main role. In Figure 19(c), it is seen that 

is not much affected by the nonuniform flow effect at low lift 
coefficients. At high lift coefficients, l4^| is slightly 
reduced by the jet flow. 
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6. Conclusions 


A lifting-surface theory for the subsonic compressible 
nonuniform flow has been developed. The theory not only 
accounts for different local dynamic pressures, but also the 
effect of velocity gradient. Comparison with limited known 
results show that the present theory is reasonably accurate. 
Numerical results indicate that there is a gain in lift if 
the wing is in a region with a positive velocity gradient. 

Based on the assumption that the suction analogy is 
still applicable in a nonuniform flow, results for a 70° - 
delta wing show that the vortex lift is enhanced by a jet flow. 

The present theory can be applied to any type of free 
stream profiles with variations in both spanwise and vertical 
directions . 
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APPENDIX A 


Integral equations for the 2-D Small-Disturbance 
Subsonic Nonuniform Flow 

The detailed derivation for equations (12) and (13) 
is given, below. 

A.l A General Integral Equation 

Applying Green's reciprocal formula to equation (3), it 
is obtained that 


= jj M + Jl (i In - p' Ji<'7 >J <>s' ( /)■ V 

S 

Now, on S ' (a small spherical surface). 

The direction of normal differential of S' points away from 
V', i.e., towards the interior of sphere. It follows that the 
second integral on the right-hand-side of equation (A.l) becomes 

I (f '3P 




as e-»-0, e ->■ 0, and S'-»- 0. Hence, the first term in equation 


(A.l), 



3M 


0 . 


Note that — fC p'c/S 


is the arithematic mean value of 


p' on S' and tends to p'(x,y,z) as^-*0. Therefore, 



4 7TP' 


Equation A.l now becomes 

I i ^ ^ ^ P5 + ^ m'> P55 j Jt'^5 

S 27 


(A.2) 



where 


4f>; = |^^-fpV|fp' 

Therefore, 

|4 = I ( r-j - if P') d5 did^ 

V 

where g = kfp'. 

From the divergence theorem, the above equation can be written 


as 



Similarly, let 


V' 



where g = jhp ' . 




C/A-4) 


Substituting equations (A. 3) and (A. 4) into equation (A.l), 


equation (A.l) becomes 


Pfx,3'Z> 


Sw 


+ a®||pW-f 

V' ( V' 






a\sy 
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Equation (A. 5) is the general integral equation for the 
small-disturbance subsonic nonuniform flow. 

Now, apply the boundary condition equation (equation 11) 
to each term (except the first term) on the right hand side 
of equation (A. 5) and perform the integration from n=“°° to 
q = 00 for the two-dimensional case. Let 


p,'= , g wx-3)YV^-5^* 

6, = 

Several terms in equation (A. 5) will be simplified separately 
in the following: 




A. 2 


.Sw 


Remember, n is positive outward away from the flow field, so 
that n = -Z for the upper surface. It follows that 


.a ( Ri' — ■3Z ^ Ri R? 


(A-V 


Substituting equation (A. 7) into equation (A. 6), it is 
obtained that 




CA- 8) 


The boundary condition requires that the flow be tangent 
to the camber surface in the thin airfoil theory. This condition 
can be written as 


Xoi-) 




- o< 




where z is the camber. To find the downwash w(x,z), equation (10) 
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is integrated to give 


ijj- ,-v ! d.x. 

— oo 


(A.IO) 

where the lower limit is chosen in such a way that w(-<»,Z)=0. 
Substituting equation (A. 8) into- equation (A. 10 ) , the following 
equation can be obtained 


w,rx,o; ^ . A 

\/«ro) ?)X ~ ^71 Jo 'X' $ J 


CAIO 


\Ao<’o) 
where 

A. 3 






(A-ny 


le-l X-1^ t 

r. I 

_ j). _L </t/; ( 

Jo V i/;^+ 

r 




r f /} 

I z 

p„ A C A. 13') 

- p, 

Substituting equation (A. 13) into equation (A. 12), it is 


obtained that 




(A- m 


Similarly, applying the boundary condition equations Eq. (9) 
and equation (11) to equation (A. 14) , it can be shown that 


Wj (X,o'> = 0^ 
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(%jO) 


= 0 . 


s- 


(A^ 15) 


Integrating from n=-<» to n=°° for equation (A. 15) , it is 
obtained that 




After applying equation (9) and equation (11) to equation (A. 16), 


the following integrals are needed: 


J-oo i^i- 3 


X dx, 


2^ X X 


J A ^ r 




(A- n) 
(A- IS) 


It follows that 


^ 2£c ^ 


2ir J) l~5' 




fAil) 


where 


Cpf5.5) = 




CA’^0) 


31 





(A^20 




*^^5“ if ^ ^ ^ ~ ^l" ^ 
E = Q.P^( 


Applying the divergence theorem to equation (A. 21), it is 
obtained that 


Ijjiiifi' -f ('■» -If ij,')* 


(A- 2 2) 


It follows that 


ft ~ 2F/ 


fj> 


c/3dj 


i A--23) 


From the boundary conditions (eq. (9) and eq. (11)), it is 
finally obtained that 


W^rxo) _-azc ^ _ -I ff 

77 ^ d : >) 
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Note: 



.pi ^ 

/X ^ [_ J. ?('l J f 

)-® < 5 ,* ^ 

f;< fjr_i>_cfy. -Ii_ 




M. 

3 [’' 

J-r, 

M iO 




f iah^ ^ f 

/5/Z-51" 





(A-! 8 ) 
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APPENDIX B 


Integral Equations for the 3-D Small-Disturbance 
Subsonic Nonuniform Flow 

To derive the wing boundary condition (eq. 15) , equation 
(A. 15) will be used again. Several terms in equation (A. 5) are 
simplified below. 

For simplicity, the following symbols are introduced: 






i -2 2 

t 1 6 Z 


2 .22 


I>,= f X-5) 

Yz=^V-pV^’ 

<3 = 

Yt=(/-pVu-p^ 


, rz=fi(/-p-^z 

^ < 5 ,= cx,-5)V/3Y-2--5Y 
, Yt'=(y-p’- 


Pi - - 4V/J I’ ■»> f R. ) '^5 


After applying (eq.9) and (eq. 11), the following equations 
can be obtained: 


W,(x,y,z)_ 3z'fx, ,/,z;dx, 


ii' = J_ / , + -fci- ) 

R,^ Ys R. ' 


■X X- ? /_L , A J_ 

R" 3YrR. R," Yz ’ 3 Yz‘ 
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It follows that 


t ' 

_ _^{( A Cp[ I ^ ( H- 

Vccy.2) Sn Yz 



B.2 




where 




a,.S) 

(B-i) 


Applying equation (11) to equation (B.6) , it is obtained that 




"X ^ 

-eo 


(X^ > 


y^z)dxt 


It follows that 


(B-V 


VJ,(X,y,^) _ I I . X -t, ; 

V^cy,z) 4t" I ft, J ^ ( 


(B.S) 


P3 - iVf 4 

where 

±£ _ p B(p 5 ^ 

R 

Applying . equation (11) to equation (B.9), the following 


integral is obtained; 




I 


( 6. 10) 


where 
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It follows that 


Wifx.y.z) , I f(((\ ,, , ,, I I 




4 Ms, 


Yt-f?-* 




(Bi^) 



where ^ ^ 

si- T^P^<y-y . ciy.’^) 

•®’[ R 


CB-13) 


(■ 6. 14) 


®^^3.tion (11) to ©quation (B.13) , it is obtained that 


Ws(x, y, ^) = - 


It follows that 



•iT 


dX, 


(B-15) 


B.5 


I. fffrp.l -Cf^-pr,-2-^fy-pfe-p X-? 
V«<y,2) 4T W i I t ^ ; 



YfR^ 


CB.it) 
fB. 17) 


By the divergence theorem 


j|v-A dv =/| A-)l dS ,g^ 

V Surface 
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CB. 19 ) 






Applying equation (B.18) to the first term on the right 
hand side of equation (B.19), it obtained that 

-l-uv =)jr» -fds 0 

' 5‘«.rj^ce 

Therefore, equation (B.17) becomes 

a- 4.1 J3 jjjj 

From equation (11) , it can be shown that 


( B.20) 



(S.VCO U 

I fff p( 


It follows that 



Wi (X,yyZ^ 

Va>cy,z) 




(6.21) 


(B. 22 ; 


B.6 Surface integrals in equations (14) and (15). 

Following the quasi-vortex-lattice method (ref. 9) , it 
is assumed that the velocity is constant for each spanwise 
vortex strip on the wing. Performing the integration along 
the bound element of each strip (Figure 4) , the double integrals 
can be reduced to finite sums through the mid-point 
trapezoidal rule. Let 
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X -^ = x-x,-r (Xj -X,) 
y- ^ = y-/i -T ( y= -y.) 

= ( x-^M/3yy-p%y5Z^ 

= + r B + c 

where 

A = ^ X2-x,)S/S\/2- /./ 

5 = -2 (;^x- x,;(Xz'XO y> 

c = ( x-x.M / 5 % /- y./+/^ 

(y-p^f = r a 4 Tb. 4 c, 

b, = -2ry-y,)^y-'y.) <^/= 

The second term on the right hand side of equation 15 can 
be written as 


Now, integrating equation (B. 23) along the bound element of 
each vortex strip and assuming the Mach number is constant for 
each spanwise strip, it can be shown that 


W 2 ( X, y^ Z) 

4 7T P® atS b.t -f C, 


, (' z^y^-/>)fx-x,-T^yz -X,)] . 7 

^0 ar*4 b,T4C, 8^+^) J 


(6.24) 



where 


(' Z ( /a - y, ) dt 
)o ai'*4b,T+C. 


- _^a±k ki 

|/2-/.llZl^ 2/(yj-y.)zl 2k/^./,;2|^ 


To integrate the second integral in equation (B.24), let 


(y^ )r -( y- y ) = V 

T = o V = - (y -y.) 

T - \ V = - / 

arH j>,t + c, = \7 % 

Ax^+ BT +c = (6.^6) 

It follows that 


r A - /.) 


2B'i/+ c')^ 


f , _i qi/-(X2-x,)z jy = X-X 

' 2 B' 9 f c' J i/ = /i '/ 


CB.^7; 


Using equations (B.25) and (B.27), equation (B.24) becomes 


1^'' 2|(si-y,>z| 2K/j-y,;z| ^ 




where 


2 ^ i . 


A' = (X.-X.) +yS (/z-yo 

B' = (fx,-x,/(y-y,)-(x-x,j> fXi -x.> </.-/. J 


c' = (Xz-x , 2 (x-x.MAz-x,)(y- x; c/z-x^ 

t (x-x.yc/z -/.)%yS zV/z- x/ 
a = ( Xz -X.) ( y- - r X- xo <r/z- xiy 
a = (/z-/,>\ t, = -2 ^y-x;yyz-x; 
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Finally, equation (B.28) becomes 

w,(x,y,z)=2;j +f] 


(B.^V 


The .first term of equation (14) can be written as 


ff fiz 


1 

Of 0^ 


Following the steps used above, let 


X - ^ = / - X, -T (Xi -X,J> 

X- 7 = x- X ~r r x^ -/,) 

(X- 5)% f r + c 

Hence , 


r B-3o; 


^,cx,y,z) 

^ (y^_y^ f __!iA±±^ 

4tt "" ‘ 1 ('^/l^c-bV/a+bk;^ ('4-/^c-b'*;c^ 


A= (x^-x.)ty5\Xi 

g = -2 (( X-X.) (X* /') 


casi; 


The second term on the right hand side of equation (14) 


is 


*^5 (X,/,z) 

- I (' ^«’>^^(y2~x) /, 

gr + cj)>^ 






(^Ac - s 


-j 1 

If -SNJiJ 


(4-Ac.~b )' 


CB.52) 
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APPENDIX C 


Integration Regions and Coordinate Transformation 
C.l Integration regions and coordinate transformation in 
the C direction are indicated in Figures 6 and 7. For 
simplicity, denote the integrand by F (5 ) . Then, 

Note that the infinite regions are actually approximated by 
appropriate finite regions in numerical calculations. The 
upstream integral becomes 

jl' CC2) 

^ = -Xci ( Xf aj) 0= ^ N, 

The second integral of equation (C.l) can be written as 

1^', (C.3) 

^iK ~ Xj, f ~ X^i i (C,3ft) 

= I- cose,) &K= ■ ' K,l,^ A/, 

where 

l\Jx = 2^^- ( NvJO) + hJ^(^)) 

and 

NW(1)+NW(2) = total aerodynamic panels on the wing in 
chordwise sections 

Similarly, the downstream integral is reduced to 


£ - fo F(®) ds (C.il-) 

ku 

^ = xt -f xt ( coso) ^ e= 2, A/3 
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C.2 Coordinate transformation in the spanwise direction 

The integration regions and coordinate transformation 
are illustrated in Figure 6. If F(n) represents the integrand, 
then, 

(C.5) 

Again, infinite integration regions are approximated by finite 
regions in the numerical calculation. 

The integration from -«> to the left wing tip can be 
reduced to 

^ - /f g; 

^ - /s ( '- 

*=11 /<= b 

<p ^ ^ ~^/z ~ Ys 

Note that is the number of integration points outside the 
wing in the spanwise direction. 

Integration points for the integration perfomed over the 
wing in the spanwise direction coincide with the ends point 
of the vortex strips. Therefore, 

^ ( I- C£>S ^ ) 

On the other hand, the integration from the right wing 
tip to “ can be reduced to 

s If Fcf'ydf (c-v 


42 


^ +/s (l-COS^; 

\ = % ^ f=% ^= ^+/S 

C.3 Integration in the vertical direction 
Assume the integrand is F(^). Then 

= C FT‘'5+X"F'y''3 

The integration from ^=-<» to r,= 0 is reduced to 

\ = -Za(l~COS'i') ^ i’=^ M, 

'^= ,V' = ^ J = 0 


where is the number of integration points in the z direction. 

On the other hand, the integration from i;=0 to 5 = 00 
becomes 

j” Feyd^ S fCf)di- cc.to) 

l-cos'f) , -V-= 

where is the number of integration points in the positive 
Z direction. 
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APPENDIX D 


Calculation of Wing Wake 


Silverstein and Katzoff (ref. 14 and ref. 15) obtained anpirical 
equations for calculating the wing wake. The itiaximum loss of dynamic 
pressure in the wake can be expressed as following equation. 

" 4- f 0,3 

c 

where is the wing section profile-drag coefficient, and x is the 

distance between wing trailing edge to the aerodynamic center of the 
tail, q^ represents the dynamic pressure of the uniform portion of 
the wake flow profile. The distribution of the dynamic pressure loss 
within the wake is q • 

% 


4m 


( P-2) 


The half width b of the wake is given by 

B.&S- c/ cl -t 0.15-)“ 

c 




Fran Equations (D-1) , (D-2) , a function of the Mach number profile 
of the wing wake can be shewn as the following expression. 


M, (l - 


^ + 0.3 


fS 


COS ( 




(p-4) 
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'fhe wake location can be determined by using the method of Chapter 7 
of reference 16. Assuming elliptic circulation distribution for wing 
in the spanwise direction, then the induced angle of attack is given as; 

^^” 5 ) 

where experimental data for Cj^ are to be used (ref. 16). 

The downwash angle can be obtained from figure 7-14 (ref. 16) by using 
cxj^. Finally, the downward displacement of the wake can be approximated by 

the following equation: 

^ ^ ^ j £ ' 

vdiere ZcQ.75sinoc for angle of attack without flow seperation. 

Zs:0.75sin (a- for angle of attack with flow seperation. 
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Table 1. Convergence Check of Numerical Integration for 
Joukowski Airfoil 


CASE 

NO. 

IGP 

.. .. . 

X 

a 


z 

a 


N1 

N2 

N3 

M 

a 

^b 

1 

5 

9.0 

8.0 

3.0 

3.0 

20 

10 

10 

10 

10 

2 

5 

9.0 

8.0 

3.0 

3.0 

15 

10 

10 

10 

10 

3 

5 

9.0 

8.0 

3.0 

3.0 

10 

10 

10 

10 

10 

4 

5 

9.0 

8.0 

3.0 

3.0 

5 

10 

5 

5 

5 

5 

5 

9.0 

8.0 

2.0 

2.0 

15 

10 

10 

5 

5 

6 

5 

9.0 

O 

00 

2.0 

2.0 

10 

10 

10 

5 

5 

7 

5 

12.0 

11.0 

2.0 

2.0 

20 

10 

20 

5 

5 

8 

1 

9.0 

8.0 

3.0 

3.0 

20 

10 

10 

10 

10 

9 

1 

9.0 

8.0 

3.0 

2.0 

10 

10 

10 

5 

5 



1 

2 

3 

4 

5 

6 

m 

00 

9 

10 



1.3235 

2.0064 



2.9887 







1.3310 



2.0108 

2.5702 

2.9185 

2.9902 

2.7468 



0.1091 



1.3498 




2.9942 

2 . 7^96 



0.1095 

B 

1.7472 

1.6748 



3.0578 

3.1057 

— 

2.8391 

2.2491 

1 .3509 


B 


1.3113 

2.0082 

2.5768 

2.9309 

3.0042 

2.7613 

2.1940 

1.3186 

0.1109 

6 

0.7253 

1.3302 

2.0196 

2.5845 

2.9365 

3.0083 

2.7643 

2.1960 

1.3198 

0.1127 

7 

0.8506 

1.3701 

2.0413 

2.5982 

2.9456 

3.0144 

2.7683 

2.1985 

1.3212 

0.1190 

8 

0.7294 

1.3249 

2.0075 

2.5683 

2.9174 

2.9896 

2.7465 

2.1863 

1 .3126 

0.1091 

9 

0.7363 

1.3345 

2.0228 

2.5875 

2.9392 

3.0110 

2.7664 

2.1976 

1.3207 

0.1116 


Notes. I indicates the pressure point locations given by Eq. (C.3a). 

Case 1 is judged to be the best solution. 

IGP=5:by least square method with 100 points of in equation (12) 
IGP=1 : with Z^=0.05 only in equation (12) 

Airfoil chord length=1.0 
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Figure 3. Vortex Element and Control Point Distribution 
Over the Wing 
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Figure 5. Integration Points and Region in Chordwise and 
Vertical Direction 
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Figure 6. Integration Points and Regions in the Spanwise Direction 
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Figure 7. Mach number profile for jet and wake streams 
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SYMBOL CONFIGURATION 

Ref. 5 

— — Present Method, M^=0.4 

Present Method, M^=0.2 



Wake 


M(z) 



mz) 


Figure 10. Lift Coefficient Ratio for an Elliptic Wing of AR= 10 . / ( 1 -Mq^ ) 
in Nonuniform Flow. Cl^ is the Lift Coefficient in ^ 

Uniform Flow, Reference Dynamic Pressure= 0.5p(0)V(0) 



( b ) 


Figure 11. Wake of the Wing at High Angles of Attack 
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SYMBOL CONFIGURATION 

Unif. (Attoched Flow) 

Unif. (Vortex Lift) 

Jet (Attached Flow) 

Jet with Vortex Lift 
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